#include <math.h>
void CORU(R, Rho)
{
	double sig3, ri3, R, CORU, rho;

	sig3 = SIG2*sqrt(SIG2);
	ri3 = sig3/(R*R*R);
	CORU = 2*PI*EPS4*(Rho*sig3)*(ri3*ri3*ri3/9-ri3/3);
}
